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Recoil “kicks” induced by gravitational radiation are expected in the inspiral and merger of black 
holes. Recently the numerical relativity community has begun to measure the significant kicks 
found when both unequal masses and spins are considered. Because understanding the cause and 
magnitude of each component of this kick may be complicated in inspiral simulations, we consider 
these effects in the context of a simple test problem. We study recoils from collisions of binaries 
with initially head-on trajectories, starting with the simplest case of equal masses with no spin; 
adding spin and varying the mass ratio, both separately and jointly. We find spin-induced recoils to 
be significant even in head-on configurations. Additionally, it appears that the scaling of transverse 
kicks with spins is consistent with post-Newtonian (PN) theory, even though the kick is generated 
in the nonlinear merger interaction, where PN theory should not apply. This suggests that a simple 
heuristic description might be effective in the estimation of spin-kicks. 

PACS numbers: 04.25.Dm, 04.25. Nx, 04.30.Db, 04.70.Bw, 95.30.Sf, 97.60. Lf 


I. INTRODUCTION 

The coalescence of spinning black holes in a binary 
system is expected to occur throughout the universe, 
on scales ranging from stellar black holes formed as the 
end-products of stellar evolution to supermassive black 
holes that lurk at the centers of galaxies. The final 
merger of such systems will produce an intense burst 
of gravitational radiation; if this radiation is emitted 
asymmetrically, as in the case of unequal masses and 
spins, the resulting remnant black hole will experience 
a recoil kick. The magnitude of this kick is impor- 
tant in a variety of astrophysical situations, such as 
the cosmological evolution of supermassive black holes 
[1, 2, 3, 4, 5, 6, 7, 8] and the growth and retention of 
intermediate-mass black holes in dense stellar clusters 
[9, 10, 11, 12, 13, 14, 15], and it also affects the ex- 
pected rates of black hole mergers for gravitational wave 
detectors [16]. Given the importance of recoil kicks in 
astrophysics, there have been numerous analytic studies 
of this phenomenon [17, 18, 19, 20, 21, 22, 23, 24, 25]. 
However, since nearly all of the recoil occurs in the regime 
of strong gravitational fields, numerical relativity simula- 
tions are essential to obtain accurate calculations of the 
kick. 

Recent breakthroughs in binary black hole simulations 
[26, 27, 28, 29] now allow extensive studies of equal mass 
nonspinning black hole mergers [30, 31, 32]. Calculations 
of the kick resulting from mergers of nonspinning black 
holes with unequal masses have been carried out for mass 
ratios q = rri 2 /mi in the range q = 1 to q = 0.25 [33, 
34, 35], with the maximum kick velocity estimated to be 
~ 175km/s for mass ratio q ~ 0.36 [35]. In addition, 
several simulations of mergers of equal mass, spinning 
black holes have been carried out [36, 37, 38]. 


Most recently, several papers have appeared that ad- 
dress the kicks obtained from the inspiral of spinning 
binaries of equal [39, 40] and unequal masses [41]. These 
papers find considerable kicks resulting from the addition 
of symmetry-breaking spins, leading to a total kick speed 
of ~ 440km/s. Using this figure, and extrapolating to an 
“optimal” (from post-Newtonian considerations) orbital 
configuration, Campanelli et al [41] suggest a maximum 
kick of ~ 1300km/s. These substantial predictions have 
opened the door to exciting new possibilities in the as- 
tronomy of supermassive black holes and galactic nuclei. 

Nevertheless the complexity of the binary orbital evo- 
lution may obscure details such as the direction of the 
final kick, and its dependence on mass ratio and spin. 
With this in mind, we present here a study of a simpler 
merger problem, which may be seen as an approximation 
to the final plunge to merger. We chose the head-on case 
as a model problem to isolate kick effects from the or- 
bital inspiral motion. Although lacking in astrophysical 
likelihood, our head-on investigations have the advantage 
of removing directional ambiguity in the kicks produced. 
In particular, we can readily test the leading order PN 
prediction that the spin and mass-ratio contributions to 
the kick should be orthogonal. 

In our investigation of the recoil kicks produced by the 
merger of spinning black holes in head-on collisions, we 
vary both the black hole spins and mass ratio. We find 
definite transverse kicks from the merger of equal-mass 
holes with spin. Furthermore, the total kick momentum 
imparted appears to scale roughly with the sum of the 
black hole spin parameter a\ + a 2 . For the cases inves- 
tigated, these kicks yield transverse kick velocities in the 
range ~ 15 — 30km/s. We have also seen longitudinal 
(along the line of motion of the holes) kicks, due to a 
non-unity mass ratio. These kicks are much more mod- 
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est in magnitude, with velocities of ~ 2 — 5km/s for the 
mild mass ratio chosen. 

This work is carried out using the moving puncture 
method [27, 28, 42, 43, 44, 45]. We describe our initial 
data in Sec. II and our numerical methodology in Sec. III. 
Code calibration and testing is presented in Sec. IV. We 
present our results in Sec. V and conclude with a discus- 
sion in Sec. VI. 


II. INITIAL DATA 

We set up initial data for binary black holes repre- 
sented as “punctures” [46]. The metric on the initial 
spacelike slice is written in the form gij = ip 4 6ij, where 
i,j = 1,2,3, with conformal factor ip = ipBL + u. The 
singular part of the conformal factor takes the form 
ip bl = 1 + zLa=i m A ; p/2|r — ta\, where the A th black 
hole has mass parameter (or “puncture mass”) ra^p and 
is located at coordinate position ta- The nonsingular 
function u is calculated by solving the Hamiltonian con- 
straint equation using a second-order convergent elliptic 
solver AMRMG [47]. 

The black holes start out at rest. We use the Bowen- 
York [48] form of the extrinsic curvature to incorporate 
black hole spin: 

3 

K A.ij ~ ~3~ [ e kim^A n A n A j + tkjmS a r ^A ri A i] j (1) 

r A 

where Ha is the unit vector in the direction of increasing 
ta and Sa is the spin angular momentum of the A th 
black hole. In all cases we take Sa to be aligned with the 
z axis, so that Sa = Sa& z - 

In all, we consider head-on collisions for seven differ- 
ent initial configurations, varying the masses and spins. 
The initial parameters of all the simulations are listed in 
Table I; all length scales are given in terms of a fiducial 
mass M (which coincides with the ADM mass of the sys- 
tem for equal-mass runs). For each run, the punctures 
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were initially placed on the y axis on the equatorial plane 
(z = 0) with the center of mass at the origin. The proper 
separation l between the holes, along the initial slice, is 
measured between the closest parts of the apparent hori- 
zons along the y axis. We note that l is not necessar- 
ily the smallest physical distance between the holes (in 
this spatial slice), as spin effects can twist the space-like 
geodesics off the coordinate axis. 

We use the horizon mass m to define the black hole 
mass ratio q on the initial slice, 

q = mi/m. 2 , (2) 


where m\ is the horizon mass of the lighter hole. The 
symmetric mass ratio is 


v = mi m 2 /(mi + m 2 ) 2 . (3) 

The horizon mass is derived from the apparent hori- 
zon’s irr'educible mass using the Christodoulou formula 
[49]:. 


2 2 S 2 

m = m ^ + 4^r 


(4) 


where mi rr = \] Aah! 167 t and Aah is the area of the 
apparent horizon. 

The spin vector S applied to each hole in the Bowen- 
York data prescription (1) is in fact the ADM angular 
momentum of that hole - the total angular momentum 
of the spatial slice as measured by an ADM integral at 
infinity if no other sources were present. Thus these pa- 
rameters represent a global quantity 

The more standard definition of astrophysical black- 
hole spin, however, is a local one. The Kerr solution is 
parametrized by the black hole horizon mass m, and a 
spin parameter a restricted to the range of values [0,m]. 
Then, a = \S\/m. We note, however, that it has proved 
impossible to bring Bowen- York spinning data to an a 
value of more than ~ 0.927m, significantly below the 
maximal Kerr value [50] . 


Even in the head-on case, with spin vectors orthogonal 
to the orbital plane, there is much freedom in the param- 
eters describing the initial data. In particular, we have 
chosen to scale all distances so that the horizon mass of 
the lighter hole, mi , is kept the same relative to the grid 
spacing (ensuring a common level of resolution of that 
hole’s features). Within this restriction, we have also 
striven to maintain the same proper horizon-to-horizon 
separation l between the holes. 


f 

III. METHODOLOGY 

These initial data sets were evolved using the moving 
puncture method as implemented in the Hahndol code 
[32]. We used standard BSSN evolution equations, with 
the addition of dissipation terms as in [51] and constraint- 
damping terms as in [52] in order to ensure robust stabil- 
ity. Our gauge conditions are 1 4-log lapse slicing and a 
hyperbolic Gamma-driver shift condition, as used in [42]; 
here we take rj = 2.0, and the initial lapse shape chosen 
is OLinit — iPo\ w ^ere ipo is the initial conformal factor; a 
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TABLE I: Simulation initial parameters in terms of a fiducial mass M. mi is the horizon mass of the lighter hole, q is the 
mass ratio defined by horizon masses. S A is the non-zero component of the Bowen-York angular momentum on each hole. 
a A = \S A \/m A is the approximate Kerr spin parameter of hole A. I is the proper horizon-horizon separation of the holes. 
NEa+_ and NEb+_ are large- separation and equal-spin-parameter variants of NE+_, as explained in the text. M A dm was 
calculated at a finite coordinate distance from the origin (60 M for all but run EQ +0 ). 


Run 

\mi, p /M m 2 ,p/M Si/M 2 

Sl/M 2 

mi/M Tin'll M 

<7 

v a\ 

U2 

ai/mi ai/mi 

l/M 

Madm/M 

EQoo 

0.5000 

0.5000 

0.0 

0.0 

0.514 

0.514 

1.0 

0.250 0.0 

0.0 

0.0 

0.0 

12.24 

1.0 

NEoo 

0.4909 

0.7478 

0.0 

0.0 

0.514 

0.771 

0.667 0.240 0.0 

0.0 

0.0 

0.0 

12.24 

1.238 

FQ + o 

0.3444 

0.5000 

0.2 

0.0 

0.514 

0.514 

1.0 

0.250 0.389 

0.0 

0.758 

0.0 

12.24 

1.0 

EQ + _ 

0.3444 

0.3444 

0.2 

-0.2 

0.514 

0.514 

1.0 

0.250 0.389 0.389 

0.758 

0.758 

12.4 ± 0.2 

1.0 

NE+- 

0.3436 

0.7140 

0.2 

-0.2 

0.516 

0.773 

0.668 0.240 0.388 0.259 

0.752 

0.335 

12.6 ± 0.2 

1.244 

NEa+_ 

0.3436 

0.7140 

0.2 

-0.2 

0.513 

0.764 

0.672 0.240 0.390 

0.262 

0.759 

0.342 

17.0 ± 0.2 

1.250 

NEb + _ 

0.3436 

0.5496 

0.2 

-0.4486 

0.516 

0.784 

0.658 0.239 0.388 0.572 

0.752 

0.730 

13.0 ± 0.2 

1.257 


steeper fall-off near the punctures than the 'ipQ 2 used in 
the evolutions of [32]. Time integration was carried out 
with a fourth-order Runge-Kutta algorithm, and spatial 
derivatives with fourth-order-accurate finite differencing 
stencils. For the outer boundary we employed a second- 
order- accurate Sommerfeld condition; since the majority 
of our simulations here last ~ 200M, we locate the outer 
boundaries at 256M to ensure that there won’t be any 
harmful effects from the boundaries. We used adaptive 
mesh refinement (AMR) implemented via the software 
package PARAMESH [53, 54], with fifth-order-accurate 
interpolation between refinement regions (as preservation 
of fourth-order accuracy in the bulk demands better than 
fourth-order-accurate interpolation at refinement inter- 
faces [44, 55, 56]). 

The momentum kick of the merged black hole is calcu- 
lated as the aggregated time-integral of the thrust , dP l /dt 
- a surface integral of the squared time-derivative of the 
radiative Weyl scalar 04 times the unit radial vector [57]: 



As this is quadratic in the waveform 04, the kick can 
be expected to be weak. To perform the angular inte- 
gration in (5), we use the second-order Misner algorithm 
described in [58, 59]. 

IV. CODE CALIBRATION 

We have established the convergence of the Hahndol 
code in prior publications [32, 34]. To these previously 
published results, we have added only one qualitatively 
new aspect: spin on the pre-merger holes. Thus we re- 
strict our discussion of convergence here to a study of 
the convergence properties of a sample spinning data set 
with equal masses and anti- aligned spins, EQ + _. 


Our main evolutions were performed at a maximum 
resolution of M/32 in the vicinity of the pre-merger holes. 
We emphasize here that, as we have maintained a smaller 
horizon mass of m ~ 0.5M in all runs, the grid spac- 
ing realized near the smaller hole will always be close to 
m/16. 



FIG. 1: Integrated transverse momentum kick A P x for the 
EQ + _ run at central resolutions of M/32 (solid/black) and 
M/40 (dashed/red). 


In our simulations, the Hamiltonian constraint ex- 
hibits greater than third-order convergence, and the 
momentum constraint exhibits second-order convergence 
through most of the duration of the runs, which appears 
to be due to our second-order- accurate initial data solver. 
At late times the apparent convergence of the momentum 
constraint drops to a rate between first and second order; 
we believe this is because of a high-frequency gauge pulse 
propagating outward until it is insufficiently resolved at 
the grid resolution achieved in the outer zones of our 




4 


coarser run. However, this gauge pulse does not sig- 
nificantly impact the thrust computation, as evidenced 
by the invariance with extraction radii demonstrated in 
Sec. V. 

In Fig. 1, we show the integrated momentum kick for 
our data set EQ + _, at our standard resolution (M/32 
in the most refined region), as compared with a higher 
(M/40 in the most refined region) resolution. The dom- 
inant uncertainty in our kick measurements is the ac- 
curacy of our extraction method. For the integrations 
performed here, the extraction method is second-order- 
accurate in the underlying grid resolution. Thus to esti- 
mate our error, we have assumed second-order accuracy 
in the integrated kick and Richardson-extrapolated our 
M/32 and M/40 runs. The difference between the result 
of this extrapolation and the M/32 kick is ~ 1.2 x 10“ 6 ; 
guided by this case, we take 2 x 10“ 6 as our error bar for 
all equal-mass runs. 

A. Bowen- York Radiation Pulses 

Bowen- York black hole initial data contains a non- 
negligible amount of “false radiation,” even after the 
Hamiltonian constraint has been solved. This radiation 
is unphysical in the sense that it is a reflection of the 
unphysical approximations, such as conformal flatness, 
made in generating the initial data (c.f. [60]). 

We can consider a Bowen- York black hole to be a per- 
turbation of a Kerr hole with (to leading order) the same 
mass and spin. The hole relaxes to a Kerr state by ra- 
diating the perturbation through its quasinormal modes 
(QNMs). Assuming the spin vector is taken in the z 
direction, the most slowly damped QNM will be the fun- 
damental (l = 2,771 = 2) mode; for a hole of mass m and 
spin parameter a = 0.7m, the damping (or “e-folding”) 
time associated with this mode is r e ~ 12m [61]. Judg- 
ing from this, we expect the BY pulses from a pre-merger 
hole of mass m to have dropped two orders in magnitude 
after 4.6r e ~ 55m. This translates to ~ 28 M for our 
equal-mass (EQ) data sets, and ~ 43 M for our unequal- 
mass (NE) data. 

This Bowen- York “pulse” is concentrated close to each 
hole’s horizon, and much of it falls into its parent hole; 
however, what radiation escapes to our detection spheres 
will mix with the physical radiation we are generally in- 
terested in. Nevertheless, since the pulse is released im- 
mediately, and is of short duration, it should be easily 
identifiable in the total waveform. 

To illustrate, we show in Fig. 2 the dominant ( l = 
2 ,m = 2) spin- 2- weighted spherical harmonic mode of 
the “outgoing” Weyl scalar 0 4 (extracted at R ex t = 
30M), for three different equal-mass configurations: 
zero-spin (EQ 00 ), single-spin (EQ +0 ) and double-spin 
(EQ + _). Despite the differences in spins, these three 
waveforms agree well in the later “physical” part of the 
signal, with stronger differences occurring in the non- 


physical BY pulse. In the zero-spin case, the BY pulse 
is almost negligible, but it is significant for EQ +0 and 

EQ+_. 



FIG. 2: Dominant (l = 2,m = 2) mode of the waveform 04 
for the three equal-mass configurations: EQ 00 (solid/black), 
EQ +0 (dashed/red), and EQ + _ (dot-dashed/green) (ex- 
tracted at coordinate distance R ex .t = 30 M). 


This same pulse will yield a non-zero contribution to 
the momentum kick. Again, this should be easy to iden- 
tify, as it will result in a plateau early in the run. Having 
identified a time at which the BY pulse and its “pseudo- 
kick” are finished, we remove this effect by beginning the 
integration of (5) from that point. 

The initial separation between two black holes is cho- 
sen large enough so that the initial BY data transients 
are separated from the merger event. This temporal dis- 
tinction means that we can choose some time ti nt after 
the passage of the transient as the starting time for our 
numerical integration of the momentum kick. 


V. RESULTS 

We group our simulations into three categories. First 
we consider the effect of unequal masses alone, where 
the kick experienced should be along the line connecting 
the holes. We then add spins to holes of equal horizon 
mass, to explore the dependence of kicks on pure spin 
anisotropy; this configuration should yield kicks orthog- 
onal to the line connecting the holes. Finally we consider 
the combination of both spin and unequal masses, to ex- 
plore the interrelation of the two mechanisms generating 
kicks. 
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FIG. 3: Longitudinal thrust dP y /dt (top) and kick AP y (bot- 
tom) for NEoo, at extraction radii R ex t of 30 M, 40 M, 50 M 
and 60 M, where the latter three have been time-shifted by 
10.6M, 21. 0M and 30.4M for alignment. 


A. Unequal masses without spin 

The head-on collision of two holes of unequal masses 
has been considered several times in the past, both with 
close-limit analysis (CLA) [62] and with fully numerical 
evolutions in 2D [63]. The CLA results (performed with 
Misner two-sheeted data) provide a definite prediction of 
preferential radiation of linear momentum; the numeri- 
cal results seemed to confirm this, and indicated where 
the CLA fails. The numerical results were for proper 
separations far less than ours; extrapolating their results 
to the larger separations we treat here indicates that a 
longitudinal kick of 0( lOkm/s) would be produced. 

In Fig. 3, we plot the y-component of the thrust dP /dt 
from Eqn (5), and its time integral from the NEoo run. 
Note that the x and 2 components of dP /dt are zero: 
the kick is longitudinal, i.e. along the axis of collision 
(y axis). The final kick velocity, given in Table II, is 
~ 2.3 — 2.5km/s. 


The final kick obtained from Fig. 3 is highly consis- 
tent across extraction radii, with a spread between the 
R ext — 30 M and R ext = 40M values of ~ 0.5%. We note 
however that the difference between R ext = 50M and 
R ext = 60 M is larger than that between R ex t = 40Af 
and R ext = 50M. This appears to be because the ex- 
traction surface at R ext — 60 M lies in a region of less 
grid refinement than the three closer surfaces. 



FIG. 4: Transverse thrust dP x /dt and kick A P x for the 
spinning, equal-mass cases EQ +0 (black/solid) and EQ + _ 
(red/dashed), at extraction radii R ex t of 30 M, where the 
EQ + _ data have been time-shifted by IM to compensate for 
gauge differences. 


B. Equal masses with spin 

We now turn to simulations of spinning black holes. 
We run two configurations (EQ +0 ,EQ + _), as indicated 
in Table I. As only one particle is spinning in the for- 
mer case, we should see only “spin-orbit” effects between 
that spin and the motion of the particle. In the latter 
case, with two spins, we might also see “spin-spin” ef- 
fects; however PN predictions for the latter are of higher 
PN order than for spin-orbit [64] 1 . 

In Fig. 4 we show the transverse kick A P x observed 
for both equal-mass data sets. We note that both the BY 
pulse and the final kick for the EQ + _ data are roughly 
double those for EQ +0 (see Table II) . This indicates that 
the total kick is close to a simple sum of the individual 
kicks. Thus any spin- spin contribution to the total kick 


1 For the EQ +0 run presented here, we have used slightly different 
choices in our numerical evolution: advection terms employed 
a third-order upwinding scheme, guardcell- filling was fourth- 
order- accurate, time stepping was performed with a second-order 
Crank-Nicholson method, no dissipation was used, and the initial 
lapse function shape was oti n it = V'bx,’ 35 m t ^ ie evolutions of 
[32]. For this data alone, the physical outer boundary was at 128 
M, but the outermost extraction here was at Rext = 30 M. Aside 
from a 1 M delay in the arrival of wavefronts at the extraction 
sphere, these minor differences have had no measurable effect on 
the measured kick. 
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is negligible in comparison to the spin-orbit term. 


C. Unequal masses with spin 

Having considered collisions between non-equal mass 
and non-spinning black holes and equal mass and spin- 
ning black holes, we combine both spin and mass ratio 
in a few cases, the NE + _, NEa + _, and NEb + _ runs. 

We show in the upper panel of Fig. 5 the thrust for the 
NE+_ case. We note that unlike the simpler cases before, 
the BY pulse and physical signal are not obviously sep- 
arated in the thrust. As a result, we may expect to have 
trouble determining where to begin the kick integrations. 

In principle, a larger initial separation would mean 
more time between the BY pulse - emitted at t = 0 and 
lasting ~ 43 M (see discussion in Section IV A) - and the 
bulk of the physical radiation, which only becomes large 
near merger. The main physical kick will certainly be 
different in magnitude between the two cases, and even 
the BY pulse magnitude may differ, depending on the 
dependence of the spectrum to the initial separation of 
the binary. Nevertheless, we can expect the BY pulse 
duration to be the same, as the duration is determined 
by the leading quasinormal frequencies of each hole, and 
these frequencies are fully determined by the mass and 
spin of each hole, which is a constant between the two 
cases. In Fig. 5, we show dP x jdt and A P x for NE+_, 
with a coordinate separation of 8 M, and NEa + _, with a 
coordinate separation of 12 M. 

The extra initial separation leads to a very similar 
physical thrust, delayed ~ 22 M relative to the original 
separation. This physical thrust begins at ~ 114 M for 
the detector at R ext = 60M. Integrating forward from 
any time up to ~ 20M before this will result in the same 

transverse kick for the NEa H case. To ensure that as 

much of the physical kick as possible is obtained in the 
smaller-separation case, we integrate the NE+_ thrust 
from tint = 114 M — 22M = 92 M at R ext = 60 M. 


In Figs. 6 and 7 we present the momentum kicks in 
the transverse (a) and longitudinal (y) directions for the 
NE+_ data set, comparing with the kicks seen from the 
NEoo and EQ + _ cases before. 

Taking the transverse direction first, we see from Fig. 
6 that the momentum kick from NE+_ is significantly 
smaller than that of EQ + _. The spin angular momen- 
tum S present in each case is the same, indicating that 
the kick is not a function of S. In contrast, the momen- 
tum kick observed from NEb+_ is much larger. 

The physical momentum kicks from these runs are pre- 
sented in Table II. We note that the ratio of trans- 
verse momentum kicks A P x between the cases NE+_ and 



i i i i i i i i i i i i i i i i 

50 100 150 200 

t/M 


FIG. 5: Transverse thrust dP x /dt (top) and kick A P x 

(bottom) for unequal-mass-with-spin data, at coordinate 
separations 8 M (NE+_; black/solid) and 12M (NEa+_; 
red/dashed). Extracted at R ext = 60M. 


NEb + _ is roughly 2/3, consistent with the PN-derived 
ratio (8). 

In addition to finite differencing inaccuracy, for these 
cases we find that the high slope in the thrust for these 
cases yields a non-negligible error associated with the 
choice in integration starting point. This uncertainty is 
around 5.5 x 10 -6 in each of A P x and A P y , and we con- 
servatively assume the full uncertainty per component to 
be 8.0 x 10 “ 6 , equivalent to ~ 2km/s. 

With this in mind, we turn to the longitudinal direc- 
tion. Here we see that the initial BY burst is comparable 
in magnitude to the physical later part. Using for consis- 
tency the same integration starting time as in the trans- 
verse direction, we present the longitudinal kicks in Table 
II. These are asterisked to emphasize our uncertainty in 
the kick magnitude. To this uncertainty, all longitudinal 
kicks are consistent with the zero-spin case, NEoo- 


D. Summary of Results 


We draw together in Table II the kicks observed in each 
of our simulations. We can develop some expectations for 
the comparative results by referencing post-Newtonian 




7 



FIG. 6. Transverse thrust dP /dt (top) and kick AP (hot- FIG. 7: Longitudinal thrust dP y /dt (top) and kick A P y (bot- 
tom) for NE 00 (black), EQ + _ (red), and NE+_ (green), and tom ) for NE 00 (black), EQ + _ (red), and NE+_ (green), and 

NEb+_ (blue). Based on integrals at R ext = 60 M. NEb+_ (blue). Based on integrals at R ext = 60 M. 


theory estimates of the radiative linear momentum loss 
due to unequal masses and spin-orbit coupling. For ex- 
ample, adapting Eqns (3.31a-b) of [64] to radial infall 
along the y axis, we find: 


Pn 

Pso 


16 dm 2 ^rriT^ 4 . 


105 mT V r J 
16 u 2 r 2 / 77i j 1 \ 4 


. 9 2mr 
—r z + | e. 


15 


^ J (a i+»2)e x , 


y> (7) 

( 8 ) 


where mr = mi +m 2 , 6m = m\ — m 2 , and r is the spatial 
separation of the two particles, with r < 0. 

In studies of spin such as the one carried out here, 
whether physical effects scale with angular momentum 
5, the Kerr parameter a = S/m or the dimensionless 
number a/m is an open question. The Post-Newtonian 
result of (8) predicts a thrust (and hence kick) that scales 
with a for each hole. We have tried to address this un- 
certainty in our choice of data sets. 

It is informative to compare the transverse kick results 
from each of our spinning data sets with each other, in 
light of post- Newtonian predictions. Looking at Eqn. 
(8), we assume that the RHS is more or less insensi- 
tive to the details of the rate of infall, the important 
length scale then being the total mass m? of the bi- 
nary. With this in mind, the remaining freedom lies in 
the spin scaling (ai + a^/mr- Relative to the EQ + _ 
case, this scaling factor is {1/2, 2/3, 2/3,1} for cases 
{ EQ +0 , NE + _ , NEa + _ , NEb+ _ } . 

In Fig. 8 we combine the post-BY-pulse part of the 
transverse thrust dP x jdt of each case, with times rescaled 



FIG. 8: Transverse thrust dP x /dt for all spinning data sets: 
EQ +0 , EQ + _, NE+_, NEa+_, and NEb+_. Time has been 
rescaled by the ADM mass for each data set, and translated 
relative to EQ + _ to line up peaks. Amplitudes have been 
scaled relative to EQ + _ according to PN predictions. 


according to the ADM mass from Table I, and translated 
so that the peaks coincide. We have also rescaled the 
thrust amplitude by the factors appropriate to the PN 
above Surprisingly, after this rescaling, the five thrusts 
fit very well, with a deviation of overall amplitude of 
~ 20%. This appears to indicate that PN predictions of 
recoil kicks have validity further into the merger regime 
that would be expected. 
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TABLE II: Final integrated momentum kicks and correspond- 
ing kick velocities. In each case, we have removed “BY pulse” 
effects through deferring integration until after the passage of 
the pulse. A P x and A P v are accurate to ~ 2 x 10 -6 for cases 
NEoo, EQ +0 , EQ + _, corresponding to errors of ~ 0.6km/s 
for equal-mass cases, and ~ 0.5km/s in unequal-mass cases. 
Cases NE + _, NEa+_, and NEb + _ suffer from additional non- 
negligible error in the choice of integration starting point. We 
conservatively estimate the total error to be ^ 8 x 10 -6 for 
these cases, translating to ~ 2.0km/s in the kick velocity. 
Asterisks mark kick values that are dominated by this uncer- 
tainty. 


Run 

tint. 

AP* 

A P y 

Vx 

Vy 

V 


(M) (10 5 ) (10 5 ) (km/s) (km/s) (km/s) 

NEoo 30 M 

51.0 

0.0 

1.2 

0.0 

2.8 

2.8 

40 M 

61.7 

0.0 

1.2 

0.0 

2.8 

2.8 

60 M 

82.9 

0.0 

1.2 

0.0 

2.9 

2.9 

EQ +0 20M 

46.1 

-5.2 

0.0 

-15.4 

0.0 

15.4 

30M 

56.6 

-5.5 

0.0 

-16.6 

0.0 

16.6 

EQ + _ 30 M 

59.6 

-10.7 

0.0 

-32.1 

0.0 

32.1 

40 M 

70.0 

-10.5 

0.0 

-31.6 

0.0 

31.6 

60 M 

90.5 

-10.8 

0.0 

-32.4 

0.0 

32.4 

NE+- 40M 

72.0 

-8.5 

1.4* 

-20.5 

3.4* 

20.8 

60M 

92.0 

-8.5 

1.5* 

-20.5 

3.6* 

20.8 

NEa+_ 40M 

94.0 

-7.9 

1.0* 

-18.9 

2.4* 

19.1 

60 M 

114.0 

-7.9 

1.0* 

-18.9 

2.4* 

19.1 

NEb+_ 40M 

72.0 

-12.3 

2.0* 

-29.3 

4.8* 

29.7 

60 M 

92.0 

-12.4 

2.0* 

-29.6 

4.8* 

30.0 


VI. DISCUSSION 

In this paper, we have addressed the roles of unequal 
masses and spins in producing recoil kicks in head-on col- 
lisions. We have chosen sample spins and mass ratios to 
explore how important spin effects are for kicks, and how 
the effects of spin and mass ratio are related in generat- 
ing kicks. We have carried out kick extractions from data 
sets with mass ratios of 1 and 2/3, and with spins on one 
or both holes. We have observed that head-on collisions 
of spinning black holes can produce significant kicks. For 
the anti-aligned spin cases we have studied, the spin-kicks 
are transverse to the direction of initial separation. The 
magnitude of the kicks, and even the thrust curves (Fig. 
8), scale with the sum of the individual black hole spin 


parameters a\ + < 22 . These kicks easily exceed the longi- 
tudinal kick produced in the case of unequal masses. To 
the accuracy available in our simulations, we find that 
the kicks due to the mass ratio and spins are indeed or- 
thogonal and independent. Both the spin scaling and the 
orthogonality of spin-kicks and mass-ratio-kicks are con- 
sistent with leading-order PN predictions for these effects 
(7,8). 

It is remarkable that the PN predictions seem to de- 
scribe our results. The kick-producing radiation derived 
in our simulations is generated in the systems’ non-linear 
mergers and the ringdowns where the assumptions be- 
hind the PN approximations clearly do not apply. In 
general terms, the PN analysis (see [64]) provides that the 
unequal-mass-kick is produced by a coupling of the mass- 
quadrupole and mass-octupole moments, while the spin- 
kick is produced by a coupling of the mass-quadrupole 
and current-quadrupole moments. At least heuristically, 
we can consider these moments even in the nonlinear 
problem. In our problem, the presence of spin does 
not seem to have a large effect on the mass moments, 
as is suggested by the spin-independence of the mass- 
quadrupole-dominated (l = 2,m = 2) waveforms, seen 
in Fig. 2. The leading effect of putting spin on the black 
holes is to scale up the current-quadrupole moment. For 
general spins, this term scales with (Si /mi — § 2 / m 2 ). 

In the more important case of inspiraling black holes, 
our results suggest that the current-quadrupole effects 
should also provide strong kicks, scaling in a similar way 
with spin. Because of the scaling with S/m we would 
generally expect the strongest spin-kicks for nearly equal 
mass mergers. Our results support the use of PN expres- 
sions (as in [39, 41]) to predict the scaling of spin-kicks for 
black hole configurations that have not yet been studied. 


Acknowledgments 

This work was supported in part by NASA grant 05- 
BEFS-05-0044. The simulations were carried out using 
Project Columbia at NASA Ames Research Center and 
at the NASA Center for Computational Sciences at God- 
dard Space Flight Center. D.C. was supported in part by 
The Korea Research Foundation and The Korean Federa- 
tion of Science and Technology Societies Grant funded by 
Korea Government (MOEHRD, Basic Research Promo- 
tion Fund). B.K. was supported by the Research Associ- 
ateship Programs Office of the National Research Council 
and the NASA Postdoctoral Program at the Oak Ridge 
Associated Universities. 


[1] D. Merritt, M. Milosavljevic, M. Favata, S. A. Hughes, [2] M. Boylan-Kolchin, C.-P. Ma, and E. Quataert. Astro- 
and D. E. Holz, Astrophys. J. 607 , L9 (2004), astro- phys. J. 613 , L37 (2004), astro-ph/0407488. 

ph/0402057. [3] Z. Haiman, Astrophys. J. 613 , 36 (2004), astro- 



9 


ph/0404196. 

[4] P. Madau and E. Quataert, Astrophys. J. 606, L17 
(2004), astro-ph/0403295. 

[5] J. Yoo and J. Miralda-Escude, Astrophys. J. 614, L25 
(2004), astro-ph/0406217. 

[6] M. Volonteri and R. Perna, Mon. Not. Roy. Astron. Soc. 
358, 913 (2005), astro-ph/0501345. 

[7] N. I. Libeskind, S. Cole, C. S. Frenk, and J. C. Helly, 
Mon. Not. Roy. Astron. Soc. 368, 1381 (2006), astro- 
ph/0512073. 

[8] M. Micic, T. Abel, and S. Sigurdsson (2006), astro- 
ph/0609443. 

[9] M. C. Miller and D. P. Hamilton, Mon. Not. Roy. Astron. 
Soc. 330, 232 (2002), astro-ph/0106188. 

[10] M. C. Miller and D. P. Hamilton, Astrophys. J. 576, 894 

(2002). 

[11] H. Mouri and Y. Taniguchi, Astrophys. J. 580, 844 
(2002), astro-ph/0208053. 

[12] M. C. Miller and E. J. M. Colbert, Int. J. Mod. Phys. 
D13, 1 (2004), astro-ph/0308402. 

[13] K. Gultekin, M. C. Miller, and D. P. Hamilton, Astro- 
phys. J. 616, 221 (2004), astro-ph/0402532. 

[14] K. Gultekin, M. Coleman Miller, and D. P. Hamilton, 
Astrophys. J. 640, 156 (2006), astro-ph/0509885. 

[15] R. M. O’Leary, F. A. Rasio, J. M. Fregeau, N. Ivanova, 
and R. O’Shaughnessy, Astrophys. J. 637, 937 (2006), 
astro-ph/0508224. 

[16] A. Sesana, F. Haardt, P. Madau, and M. Volonteri, As- 
trophys. J. 623, 23 (2005), astro-ph/0409255. 

[17] A. Peres, Phys. Rev. 128, 2471 (1962). 

[18] J. D. Bekenstein, Phys. Rev. D 183, 657 (1973). 

[19] M. J. Fitchett, Mon. Not. R. Astron. Soc. 203, 1049 
(1983). 

[20] M. J. Fitchett and S. Detweiler, Mon. Not. R. Astron. 
Soc. 211, 933 (1984). 

[21] I. Redmount and M. J. Rees, Commun. Astrophys. 14, 
165 (1989). 

[22] A. G. Wiseman, Phys. Rev. D 46, 1517 (1992). 

[23] M. Favata, S. A. Hughes, and D. E. Holz, Astrophys. J. 
607, L5 (2004), astro-ph/0402056. 

[24] L. Blanchet, M. S. S. Qusailah, and C. M. Will, Astro- 
phys. J. 635, 508 (2005), astro-ph/0507692. 

[25] T. Damour and A. Gopakumar, Phys. Rev. D 73, 124006 
(2006), gr-qc/0602117. 

[26] F. Pretorius, Phys. Rev. Lett. 95, 121101 (2005), gr- 
qc/0507014. 

[27] M. Campanelli, C. O. Lousto, P. Marronetti, and 
Y. Zlochower, Phys. Rev. Lett. 96, 111101 (2006), gr- 
qc/0511048. 

[28] J. G. Baker, J. Centrella, D.-I. Choi, M. Koppitz, and 
J. van Meter, Phys. Rev. Lett. 96, 111102 (2006), gr- 
qc/0511103. 

[29] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 
Rev. D 73, 061501 (2006), gr-qc/0601091. 

[30] J. G. Baker, J. Centrella, D.-I. Choi, M. Koppitz, and 
J. van Meter, Phys. Rev. D 73, 104002 (2006), gr- 
qc/0602026. 

[31] A. Buonanno, G. B. Cook, and F. Pretorius, unpublished 
(2006), gr-qc/0610122. 

[32] J. G. Baker, S. T. McWilliams, J. R. van Meter, J. Cen- 
trella, D.-I. Choi, B. J. Kelly, and M. Koppitz (2006). 

[33] F. Herrmann, D. Shoemaker, and P. Laguna (2006), gr- 


qc/0601026. 

[34] J. G. Baker, J. Centrella, D.-I. Choi, M. Koppitz, J. van 
Meter, and M. C. Miller, Astrophys. J. Lett. 653, L93 
(2006), astro-ph/0603204. 

[35] J. A. Gonzalez, U. Sperhake, B. Briigmann, M. Hannam, 
and S. Husa, unpublished (2006), gr-qc/0610154. 

[36] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 
Rev. D 74, 041501 (2006), gr-qc/0604012. 

[37] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 
Rev. D 74, 084023 (2006), astro-ph/0608275. 

[38] M. Campanelli, C. O. Lousto, Y. Zlochower, B. Krishnan, 
and D. Merritt (2006), gr-qc/0612076. 

[39] F. Herrmann, I. Hinder, D. Shoemaker, P. Laguna, and 
R. A. Matzner (2007), gr-qc/0701143. 

[40] M. Koppitz, D. Pollney, C. Reisswig, L. Rezzolla, 
J. thornburg, P. Diener, and E. Schnetter (2007), gr- 
qc/0701163. 

[41] M. Campanelli, C. O. Lousto, Y. Zlochower, and D. Mer- 
ritt (2007), gr-qc/0701164. 

[42] J. R. van Meter, J. G. Baker, M. Koppitz, and D.-I. Choi, 
Phys. Rev. D 73, 124011 (2006), gr-qc/0605030. 

[43] M. Hannam, S. Husa, D. Pollney, B. Briigmann, and 
N. 6 Murchadha (2006), gr-qc/0606099. 

[44] B. Briigmann, J. A. Gonzalez, M. Hannam, S. Husa, 
U. Sperhake, and W. Tichy (2006), gr-qc/0610128. 

[45] U. Sperhake (2006), gr-qc/0606079. 

[46] S. Brandt and B. Briigmann, Phys. Rev. Lett. 78, 3606 
(1997), gr-qc/9703066. 

[47] J. D. Brown and L. L. Lowe, J. Comput. Phys. 209, 582 
(2005), gr-qc/0411112. 

[48] J. Bowen and J. W. York, Phys. Rev. D 21, 2047 (1980). 

[49] D. Christodoulou, Phys. Rev. Lett. 25, 1596 (1970). 

[50] S. Dain, C. O. Lousto, and R. Takahashi, Phys. Rev. D 
65, 104038 (2002), gr-qc/0201062. 

[51] P. Hiibner, Class. Quantum Grav. 16, 2823 (1999). 

[52] M. D. Duez, S. L. Shapiro, and H.-J. Yo, Phys. Rev. D 
69, 104016 (2004), gr-qc/0401076. 

[53] P. MacNeice, K. Olson, C. Mobarry, R. de Fainchtein, 
and C. Packer, Computer Physics Comm. 126, 330 
( 2000 ). 

[54] http : //ct . gsf c . nasa . gov/paramesh/Users_manual/ 
amr.html. 

[55] B. Imbiriba et al., Phys. Rev. D 70, 124025 (2004), gr- 
qc/0403048. 

[56] L. Lehner, S. L. Liebling, and O. Reula, Class. Quantum 
Grav. 23, S421 (2006), gr-qc/0510111. 

[57] E. T. Newman and K. P. Tod, in General Relativity and 
Gravitation: One hundred years after the birth of Albert 
Einstein, vol. 2, edited by A. Held (Plenum Press, NY, 
1981), pp. 1-36. 

[58] C. W. Misner, Class. Quantum Grav. 21, S243 (2004), 
gr-qc/99 10044. 

[59] D. R. Fiske, J. G. Baker, J. R. van Meter, D.-I. Choi, 
and J. M. Centrella, Phys. Rev. D 71, 104036 (2005), 
gr-qc/0503100. 

[60] M. Hannam, S. Husa, B. Briigmann, J. A. Gonzalez, and 
U. Sperhake (2006), gr-qc/0612001. 

[61] E. Berti, V. Cardoso, and C. M. Will, Phys. Rev. D 73, 
064030 (2006), gr-qc/0512160. 

[62] Z. Andrade and R. H. Price, Phys. Rev. D 56, 6336 
(1997), gr-qc/9611022. 

[63] P. Anninos and S. Brandt, Phys. Rev. Lett. 81, 508 



(1998), gr-qc/9806031. 

[64] L. Kidder, Phys. Rev. D 52, 821 (1995), gr-qc/9506022. 



